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1—1 Abstract 

> 

We study the collinear splitting functions needed for next-to-next-to-leading order 
calculations of jet production in the vacuum and in dense QCD matter. These splitting 
functions describe the probability of a parton to evolve into three-parton final state and 
■^J- are generalizations of the traditional DGLAP splitting kernels to a higher perturbative 

order. Of particular interest are the angular distributions of such splitting functions, 
which can elucidate the significance of multiple parton branching for jet observables and 
guide the construction of parton shower Monte Carlo generators. We find that to 0(ot 2 s ) 
both the vacuum and the in-medium collinear splitting functions are neither angular 
ordered nor anti-angular ordered. Specifically, in dense QCD matter they retain the 
characteristic broad angular distribution already found in the 0{a s ) result. 



1 Introduction 

To understand the various aspects of an experimental measurement, the LHC program relies 
on simulations from parton shower Monte Carlo codes like PYTHIA pQ and HERWIG [2]. 
One of the ingredients of such event generators are the collinear Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) splitting functions [31 IH E]. These functions describe to leading 
order 0(a s ) collinear radiation on a distance scales of the order of the inverse transverse size 
of a typical jet. Coherent branching effects J6J [7j are encoded in some of these generators, like 
HERWIG, in a way that forces subsequent branchings to happen at smaller angles than the 
previous ones. This property is known as angular ordering. Typical parton showers resum 
leading order collinear large logarithms, while the angular-ordered parton shower have been 
argued to include leading order infrared logarithms. 

While implementations of angular ordered parton showers have been phenomenologically 
successful, they suffer from a conceptual inconsistency. Angular ordering is directly applied to 
the collinear splitting functions. The coherent branching is correct in the soft gluon emission 
limit from a hard iV-parton final state. Conversely, the collinear (DGLAP) splitting functions 
are derived in the small angle approximation, factorize from the hard scattering and have no 
knowledge of the global event structure. In this paper we are interested in characteristics of 
well-separated and energetic jets. In order to understand the properties of highly-collimated 
showers, and in particular the questions of angular ordering, angular anti-ordering, or lack 
of any ordering, we investigate the higher order 0(a 2 s ) collinear splitting functions. In the 
vacuum, all such 1 — > 3 parton branchings have already been calculated [HI E] • We use Soft 
Collinear Effective Theory (SCET) [THl EH CGI H3] to demonstrate, on the example of the q — > 
ggq splitting, that it recovers these results. Our main focus, however, are the medium-induced 
splitting processes, where the characteristic large-angle radiation pattern, first understood in 
the soft gluon approximation limit [14J, is the theoretical basis for interpreting jet production 
in heavy ion collisions at the LHC. Naturally, the question of what effects, if any, multiple 
branchings may have on the angular distribution of an in-medium shower is a very important 
one. So far, possible qualitative features of the gluon bremsstrahlung have been discussed on 
the example of a dipole antenna model [T5l IT6] . Actual calculations of 1 — > 3 parton branching 
in dense QCD matter are absent in the literature. We use Soft Collinear Effective Theory with 
Glauber gluons PH HEJ E3 HH [2D 1221 123] (SCET G ) to derive the splitting function of a quark 
to emit two gluons and study the angular distributions in such higher order medium-induced 
splitting. This information might lead to important insights as to what corrections arise to 
fixed order and/or resummed calculations [2H EH EH1 EZl EHl [291 EQl EH [32] and Monte Carlo 
simulations [331 EH ESI EHl E7] of jet observables in heavy ion reactions. It will also help 
interpret the exciting experimental measurements with jet final states in heavy ion collisions 
at the LHC, see for example [331 ESI EH SU ■ 

Our paper is organized as follows. In section [2] we derive the q — > ggq splitting function 
using SCET. In section [3] we review the basics of Soft Collinear Effective Theory with Glauber 
gluons. In section[4]use SCETq to derive the first order in opacity splitting function of q — > ggq. 
In section [5] we study the angular distributions of vacuum and medium-induced 1 — > 3 splitting 
functions. We conclude our paper in section [6j 



2 The splitting function q — > ggq in the vacuum 

In this section we calculate the vacuum splitting function for q — > ggq\j We use SCET 
and show, that our SCET calculation yields the same result as obtained in the collinear 
approximation of massless QCD, Ref. [8]. 

We follow the method used in Ref. [211] for radiative energy loss. We use the light-cone 
gauge which allows us to work directly with physical transverse polarization vectors. For aThe 
gluon polarization vector e in light-cone components is 



\n-e,n-e } e±\ 



n-p 



(1) 



where we have used the gauge conditions e(p)-p = 0,n-e = 0. The spin sum matrix of these 
transverse polarization vectors is 

£ e\e{ = $*. (2) 

gluon polarizations 

i and i' are indices in Euclidean three-space. 

2.1 Notation and kinematics 

We consider a general hard scattering amplitude J that creates, apart from other partons, the 
collinear parent quark with momentum p$. This parton subsequently emits two gluons with 
momenta p\ and pi-, and the final-state quark emerges with momentum p 3 . The Feynman 
diagrams for matrix elements 2 Ain and <A4„ i. 2 in SCET are shown in Figure 1 . These matrix 
elements equal to 

Mf = Xn, P0 J, (3) 

M { °l 2 = g 2 e\\el\Xn, P3 r^J, (4) 

where Xn is the gauge invariant collinear quark field of SCET and the expression for T e s follows 
directly from the Feynman rules of SCET. The squared matrix elements can be written as 



spin, color 



£ |M 0) l a = Tr (**'*> J(po)J(po)) , (5) 

in, color ' 

£ \Mnl 2 =£ 4 Tr I -n-p 3 J(po)J(p ) Po) , where (6) 

spin, color ' 

po= E VCrSftVrS?. (7) 

The trace in the equations above is over spin and color indices. To obtain Eq. (fsl) we have 
used Eq. rt2J) to sum over the gluon polarizations. If 

1 Note that to O(o^) the kinematic variables, such as the light cone momentum fractions and transverse 
momenta do not factorize. 

Subscripts ^ refer to the fact that we deal with vacuum splittings in this subsection. 
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Figure 1: Definition of the matrix elements M„ and M^ 2 and the Feynman diagrams in 
SCET that contribute to the splitting q — >■ ggq. 

PO = tp (I D irac) (Icolor) , (8) 

with ip a real numbeiQ the n + 2 parton matrix element and the n parton matrix element 
can be related 

EK | 2 | 2 = £(^we|M 0) | 2 . o) 

The sums are over spin and color. The splitting function^ for q — > ggq 

tfU*) = ^f^o, (10) 

where Zi = n-pi/n-(pi + p 2 + Pz) — n-pi/n-p , and Si 23 = (pi + P2 + P3) 2 , is calculated by 
substituting Eq. rt8J) into Eq. ^ and comparing with Eq. ^ and Eq. ([9]). The factorization 



formula Eq. d9J) and the splitting function Eq. ( 10 ) were first derived in Ref . [8] . Note that if 
not for Eq. (|9]), a process independent splitting function could not be defined. 
In section 12.21 we use the vectors 

where the four- vectors Qi and Q 2 are linear combinations of Pi,P2,Ps- They are related to 

Sij = (Pi +Pj) 2 via 

s i3 = ZiZ 3 U piP3 , s 2 3 = z 2 z 3 U p2P3 , s 12 = ziz 2 U pip2 . (12) 

These relations are needed to compare our results to Ref. [H] • Note that out of the six transverse 



vectors that appear in vacuum Feynman diagrams (see section 2.2), U PljP3 , U P2iP3 , U Pl]P2 , 
U P2>Pl+P3 , U piiP2+P3 , U pi+P2iP3 , only two are linearly independent PI They all can be written as 
linear combinations of U PliP3 , U P2iP3 and, moreover, any product of these six vectors can be 
written as combination of s i3 , s 23 , s 12 with coefficients that depend on zi, z 2 , z 3 . 



3 This is the case if the jet has been created by a pure QCD interaction. 

4 The brackets in this notation indicate that we sum over the initial quark polarization. In all splittings 
that we consider in this paper we do the same. Spin correlations in the vacuum splittings have been studied 
inRefs. [HI ID] - 

5 Due to boost invariance along the collinear direction only two of the three momenta Pi,P2,P3 are inde- 
pendent. 



2.2 Individual contributions of diagrams 

In this subsection we calculate the diagrams that contribute to the effective vertex r e fj, shown 
in Figure [TJ We use the SCET Feynman rules, the polarization vectors in the light-cone gauge 
given in Eq. pj, and the relation 

f ±7 { = -$* - ^' 3 £ 3 , where i,j = 1, 2, and S 3 = ( ° ° 3 J . (13) 

The explicit form of the matrix S 3 is only valid in the Weyl representation. Its properties, 



C,= 






Sl23 


Co = 


1 






S123S13 


c 3 = 


1 






■S123S23 


n. — 


1 



£ 3 t = X 3 and (X 3 ) 2 = 1, which we use in section 2.3, however, are independent of the 
representation. For the contributions from each of the four diagrams to T e g we obtain Cj T* 1 * 2 , 
with 

T\ 112 = [T ai T a2 O^ 2 + T a2 T ai Oif ) , (14) 

p«i«2 _ \ r rpa\rpa,2 t jjr'i t tJ2 r)i\iihfa (TK\ 

h ,h 

p«l«2 _ \ ^ rpa 2 rpa 1 TTjl T tJ2 /")hhjlJ2 ( ^ R\ 

L 3 — 2-^t lJ P2,P3 yJ Pl,P2+P3 LJ 3 ) l i0 J 

ii j'2 

p*i«2 _ \ ^[T^ai r T a 2l Tpi TT J2 /^\hi2JiJ2 (-\7\ 

i 4 — /J J ) J J U pi,p 2 U Pi+P2,P3 W 4 • \ 1 ') 

Sl23 s 12 ^T 

JU2 

The operators Oj depend only on z±, z%, z 3 and £3 and are defined as 

D — ( r\i\i.2. nii!2. /^)iii2JiJ2. r^-inhh . r^i^hh \ 
\ W ^a ) w lb ' ^2 ) ^3 > W A J 1 

Q = ( S ili2 ' £ ili23 - $*ih$i2J2. §ilJ2fii2Jl. fihi2fijlJ2. fihjl £ i2J23. fii2J2 £ hjl3. fiill2 £ jlJ23 j 

Oi = J2Qj M Ji- ( 18 ) 

3 

The 8x5 matrix M is equal to 

Ma = (19) 

/ 1 _ 2(zi-z 2 ) 1 , 2(zi-z 2 ) n n \ 

/ i-* a (z 1+ z 2 y i- Zl ■+■ ( Zl + Z2 )2 u u u \ 

— g 5 - -£?- 

-1 + 22 -1 + 21 

(2 - z 2 )( Zl + 2z 3 ) z x z 2 -2z 2 (l + z 3 ) 

Zl z 2 (2 - zi)(z 2 + 2z 3 ) -2zi(l + z 3 ) 

-2iz 2 —Z1Z2 



2ziz 2 (l+z 3 ) 
zi+z 2 



-iz 2 (zi + 2z 3 )£ 3 i(-2 + Z!)z 2 S3 2iz 2 (zi + z 2 )£ 3 

izi(-2 + z 2 )S 3 -izi(z2 + 2z 3 )E 3 2izi(zi + z 2 )£ 3 

-2i((-l + z 1 )z 2 + 2iz 3 )E 3 2iz 2 E 3 / 



j 1 



Note that Eq. (14) and Eq. (17) do not coincide exactly with the corresponding diagrams 1 



and 4 in Figure [T] since we have rearranged terms in these equations. The terms proportional 
to 2(zi — Z2)/(z\ + z 2 ) 2 in the first row of Mji are included into operators 0\ a and Ou and 
therefore contribute to diagram 1. In Figure [TJ these terms are contained in diagram 4. This 
rearrangement does not effect the sum of diagrams 1 and 4 and is especially convenient for 
the in-midium calculation. 



2.3 Result for the vacuum splitting 

As one can see from the previous subsection, the sum of all four diagrams can be written in 
the following form, 



r^ 2 = Y^ CfclT 2 = [T ai T a2 (a\)[ 2 + ia 2 )\ 2 S 3 ) + T a2 T ai (a\^ 2 + ia^ 2 S 3 )] . 



fe=i 



E e r «f + <r ^) 



(20) 



,(o) 



In general, e!- is a list of color operators and c^ 1 ] 2 , c^ 1 ! 2 are complex numbers. In our case 
e\ = T ai T a2 , e 2 = T a2 T ai and aij, a 2 j are real functions which can be extracted from 
Eq. (ph-Eq. (p|. 



Squaring r eff yields 



A) 



E 7 (rS 9 ) t 7°rS 8 = E< e ? > l e i 0) > ( Re «r«u + <?"***)) 



*1.*2=1 



/J 



— S 3 Im (q;*j,-q;2,j — <^2,f a i,j)) ■ 



(21) 



In the equations above, the dot between a's is a shorthand notation for summing over ii,i 2 , 
for example a\ .■/•aij = ^ j =i a^ 1 * 2 *^ 1 * 2 . Note that The Gram matrix for the two basis 

color operators e^ 



(0) 



(eflef) 



c 



c F (C F - Of) 



c F (C F - Si) c 



i2 



X I, 



color j 



(22) 



is symmetricn and proportional to unity in color space. Since the tensors a are real the 



imaginary part in Eq. (21) vanishes and the Dirac part is also proportional to unity. As an 



immediate result we obtain that Eq. rt8|) holds, 

po = ^o Ofoirac) (Icoior) , where (23) 

4>o = C F [Cp ((«!,! + a 2 ,i) 2 + («i, 2 + a2,2) 2 ) - C A (ai,i-a 2 ,i + a 1>2 -a 2 , 2 )] . (24) 



,(°)u(°) 



We used that the Gram matrix (e^-, |ej- ) is a symmetric matrix to arrive at Eq. (21 ) 



As before, it is understood that a square or a product contains summation over ii,i 2 . 



Substituting Eq. (24) into Eq. (10) yields 



(P 9m ) = Cl{P^l q3 ) + C F C A {P^l A ), where 



/p(nab) \ 



^3S 123 



-ai,r«2,i - ai,2-a2,2j 



(25) 
(26) 

(27) 



919293' i 

In the abelian part, the contribution from I\ cancels exactly in the sums a\ t \ + «2,i and 



ai,2 + «2,2 because of the color commutator in Eq. (17). The remaining three diagrams result 
in 



/ p(ab) \ 

V 0109 53/ 



'123 



919293' 9e e r r 

The non-abelian part reduces to 



*&(! + 2 i) , s i23 23(1 - *i) + (1 - 22) 



«13 



Z1Z 2 



■§23 
Sl3 



+ (1^2). 



(28) 



/p(nab) \ 
\ 919203/ 



[2(zis 23 - z 2 si 3 ) + (zi - z 2 )s 12 y 



4(zi + ztfsl 



+ T + 



'123 



12 



4 2si9S 



1 + z\ 1 + (1 - z 2 ) 



12*13 



2 2 



l-* 3 



'123 



23(1 + zf ) S123 ( Zi{2- 2z\ + z 2 ) - 2:2(6 - 6z 2 + z| 



4Si 3 S 2 3 ZiZ 2 

S123 



+ 



2,s 



12 



1 - 2; 2 ) 3 + z\ - z 2 



zs(l - Z\ t 



z 2 (l- z 3 ) 



22 



2.s 



13 



z 2 (l - z 3 ) 



Z\Z 2 



+ (!<-► 2). 



(29) 



Both Eq. (28) and Eq. (29) are in agreement with Ref. [8] for e — 0. 



2.4 Cascade approximation for the two gluon splitting function 

In this section we derive an approximation for the q — > ggq splitting function based solely 
on our knowledge of 1 — > 2 splitting functions. We refer to this approximation as " cascade" , 
since it is closely related to how parton shower generators would approximate such higher 
order collinear splitting function. We start from the definition of arbitrary 1 — > 2 and 1 — > 3 
splittings 



E K+. 

spin, color 



29 



E \M 

spin, color 



(0) 
n+2 



Sjk 

4/ 

s 2 

s jkl 



J -(pS 3k \Pi,p k }) E \ M 

spin, color 



(0)1 
n 



>(0) 



( P i4iki\Pj>Pk,Pi]) 



E 

spin, color 






(30) 
(31) 



Next we calculate an expression for 



M 



(o) 

n+2 



by iteratively applying Eq. (30) twice and sum 



ming over all possible branching sequences that produce two gluons. Comparing this expression 



to Eq. (31), we obtain the cascade 1 — > 3 splitting function 

( P q%g[P2,P3+Pl])(Pi%q[Pl,P3\) (Pi%q[Pl,P3+P2])(Pq%q[P2,P3\) 



s 123 I 



«13 -523 

\ P q%q [Pi + P2 , gg] ) (-Pff flfl [Pi ; gg] ) . ( .-,,,. 

Sl2 



In section 5.1 we compare the cascade approximation to the full splitting. The cascade formula 



omits certain interference terms, which are contained in the full 1 — > 3 splitting. 

3 Soft Collinear Effective Theory with Glauber Gluons 

When a highly energetic parton traverses dense QCD matter, the perturbative QCD approach 
can be used to describe its elastic and inelastic interactions and the formation of an in-medium 
parton shower. In this approach, the medium can be modeled as consisting of effective scat- 
tering centers that provide a color-screened Coulomb potential, which serves as a background 
field for the partons that travel through the medium |42j . Consequently, the processes that 
characterize the evolution of a parton shower in strongly-interacting matter, which can be 
cold nuclear matter or a Quark-Gluon-Plasma (QGP), can be divided into two categories. 
In the first category the familiar soft and collinear splittings appear at leading order at high 
energies, analogously to the vacuum case. The second category involves elastic scattering 
with the medium quasi-particles. The first type of processes are described by the known Soft 
Collinear Effective Theory. For example, small angle collinear radiation in the parton shower 
is correctly captured by the SCET Lagrangian. However, the elastic scattering off of medium 
quasi-particles forces us to go beyond traditional SCET. 

Soft Collinear Effective Theory with Glauber gluons (SCETq ) is an effective theory appro- 
priate for describing parton shower formation in the ambiance of dense QCD matter and the 
corresponding jet observables in heavy ion collisions. In addition to the interactions of SCET, 
it has interactions of collinear quarks [T7J [20] and collinear gluons [20] with t— channel off-shell 
gluons with momentum scaling^] (A 2 , A 2 , A), which are usually called Glauber gluons. So far, 
the soft gluons have been neglected and SCETq contains only interactions of collinear fields 
with Glauber gluons. Because this mode is off-shell, the proper description for it is to treat 
the source field and the Glauber gluon as a background field. Thus, based on the assumptions 
for the momentum scaling of the source, as well as on the gauge fixing, one can derive the 
scaling of the background field created by the source. With this scaling at hand, it is a matter 
of putting this background field into the covariant derivative of the SCET Lagrangian and 
extracting the Feynman rules of SCETq . The resulting Lagrangian of SCETq is 



£-SCETG{£,n, A n , A G ) — £sCET(£,n, A n ) (33) 



7 This statement is correct for the static and collinear sources. However for the soft source the correct mode 
is (A, A 2 , A) [H]. 



+*£ 



p,p 



e ' [£,n,p>T —in,p — if A np ,A np g uX n-p j n-A G (x 



The details of the Lagrangian depend on the type of the source and the gauge fixing condition. 
In [20] different types of such choices have been considered. The Lagrangian above corresponds 
to the static source with the momentum scaling p SOU rce = Mv + k, where the mass of the source 
particle M — > oo and k ~ (A, A, A). As for the gauge choice, it is the hybrid gauge, when the 
collinear gluons are quantized in the light-cone gauge and the potential off-shell Glauber gluons 
are quantized in the covariant gauge. This choice is simple for two reasons. First, the number 
of Feynman rules and their structure is minimal with this choice. For example, one can 



compare the two terms in Eq. (33) with similar Lagrangians of SCETq derived in [20] using 
covariant or light-cone gauge. Second, the two Wilson lines that appear in matrix elements 
of the effective theory, the collinear Wilson line and the transverse gauge link, do not add 
additional diagram^] [20] . 



The Lagrangian in Eq. (33) contains the background field in position space. The Feynman 
rules of such a Lagrangian contain the Fourier transform of the vector potential n-Ac{x) — > 
J d 4 q/(2ii) 4: e lQX v(q). For the static source v(q) = 2Tc5(q°)v(q) due to the fact that the recoil 
energy is negligible. This formally follows from time independence of the background field 
A G (x). 

Since every appearance of the Glauber gluon interaction leads to the integral over the 
Glauber gluon momentum, the following notation will be useful: 

d$i = J^ e^ v( Qi ), d$ J± = |^ e-**^ «(<Ui), (34) 

where 5xi = Xi — xo, xq is the space-time position where the jet was created, and Xi is the 
space-time position of the interaction with the medium quasi-particle i. The transverse part 
of the four- vector 5xi is defined as 5xj. The relation between these two definition is simple: 

dS^dSii^J^e*^**, (35) 

2tt 

where 5zi = 5xf. Finally, in order to relate the cross section to physical observables as elastic 
scattering length and cross sections, we use: 

^ 6 F 6 A , for qg -> qg 

In the equation above, C^-R) and C^iT) are the quadratic Casimirs of the incident parton 
and target (source) representations. d,A = 8 is the dimension of the adjoint representation. 
The formula above is valid in the high energy limit and neglecting the masses of the partons. 
As a result one can read out the value of v(q±) = 47ra s /(/i 2 + q^_). 



3 This was derived specifically for the static source. 
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Figure 2: Single Born diagrams. See text for explanation of the different topologies. 

4 The g — >• p<7<7 splitting function in dense QCD matter 

In this section we calculate the q — > ggq splitting function in the medium to first order in 
opacity, using SCET G and keeping the full z\, z 2 , z 3 dependence. First order of opacity contains 
single Born diagrams, representing interactions of the propagating system at longitudinal 
positions xf . It also contains double Born diagrams, which can be viewed as the contact limit 
x| — y xf of 2 interactions. In the first type of interactions one Glauber gluon is exchanged 
in both the matrix element and the complex conjugate of the matrix element. In the second 
type two Glauber gluons at the same point are exchanged either only in the matrix element or 
only in the complex conjugate of the matrix element. The organization of the opacity series 
is independent of the propagating system, for more details see Refs. [HH S3]- 

The calculation in medium is very similar to the one in vacuum and we use many definitions 
of section [2] in this section. In particular, since Glauber gluons do not carry large momenta, the 
entire part that depends solely on zi,z 2 , z 3 is identical in the vacuum and medium calculations. 



Thus, we use the same operators Oj given in Eq. (18) as well as the same matrix given in 
Eq. (ph. 



4.1 Single Born diagrams 

All single Born graphs are shown in Figure [2j Graphs 1-4 have the same collinear structure 
as the vacuum graph 1 in Figure [TJ and we refer to them as of topology 1 . Graphs 5-9 are of 
topology 2, graphs 10-14 are of topology 3, and graphs 15-19 are of topology 4. The amplitude 
of an arbitrary single Born graph k with 1 < k < 19 looks like 

M k x) = -g 2 4 1 4 2 Xn, P (J d® ± C k rf 2 4 1} ) J- (37) 



The minus sign in Eq. (37) cancels when squaring the matrix element. It is chosen for conve- 



nience, since it leads to a color operator matrix with more positive than negative numbers, see 



Eq. (95). In the remainder of this subsection we give detailed expressions for the longitudinal 



integrals I k , the factors C k and effective vertices Y k 112 . 
The longitudinal integrals are defined as 



/f = / ^^ A 9(Qu q) .- \(QN k , q), (3? 



where the integrand contains the product of all propagators with momentum Qi — q, that 
depend on the medium transfer momentum q 

A g (p, q) = [0(p, q ± ) - q~ + ie/n-p] , (39) 

n(p, q± )= P -- {p± ~ + q±)2 . (40) 

For diagram k the number of q dependent propagators is iV^. For single Born diagrams Nk is 
at least iV min = 1 and at most N max = 2 for topology 1 and iV max = 3 for the other topologies 
2,3,4. 

Performing the integrals yields 



" mm 



i Ji^), if N k = N n 

4 1] = I i 2 (n u n 2 ), if N k = N min + i } ( n ) 

[13(^1,^2,^3), if N k = N min + 2 
where Qi = ft(Qi, q_i_) and 

J^fii) = -ie iQl5z , (42) 

i 2 (n 1 ,n 2 ) = i - _ n — , (43) 

(pif22(5^ piQiSz piQ:jSz pjni(52:\ i 

^w^ n^r) n^W 2 < 44 > 

For details on the longitudinal integrals of single and double Born graphs, see Appendix [C] 

10 



The factors Cu and effective vertices Tl 1 * 2 are 



1 k — e k u la "+" e fc ^16 ' °' 



^, if N k = N aia 

4, if N k = N min + l 



(45) 



P*l«2 _ „ TTjfl TTJ2 ffiiinh n 

L k - e fc u p fcl ,p fc2 u p fc3 ,p fc4 u tfc ' °fc 



2 2 , if Nk — N m i n 
a -- , if A^ fc = N min + 1 

Pfe 12 n 'P0' 

if N k = N min + 2 



n-(pk 1 +Pk 2 ) n -P0' 



(46) 



where Pk 1 ,Pk 2 are the two four- vectors that come out of the second collinear splitting and 
similarly Pk 3 ,Pk 4 are those coming out of the first splitting. Since Uq 1i q 2 is antisymmetric 
under exchange of its arguments Q\ -H- Q2, we need to define the order of the arguments: 
For q — >■ gq splittings, the gluon momentum is the first argument of U followed by the quark 
momentum; for g — >■ gg splittings, the momentum containing p\ is the first argument of U 
followed by the momentum containing p 2 . 



The color operator for single Bor n amplitudes e^ is provided in Appendix D in the basis 



of six elements e\ , see section 



4.3 



For topology k = 1 there are two color operators per 
diagram, with indices ^ and w, while for the other topologies there is only one color operator 



per diagram, consistent with the notation in Eq. (45) and Eq. (46). 



Note that Eqs. (|45|) and (46) are very similar to the corresponding vacuum equations, 
Eq. (14)— Eq. ([l~7|)p| We also have used the same rearrangements between topology 1 and 4 
as in the vacuum case, mentioned at the end of section 12.2 



Even though we provided all rules necessary to evaluate any single Born graph in this 



subsection, we also summarize all values for C*., U 



Pk x ,Vk 2 1 ^Pk 3 ,Pk 4 



rW 



in Appendix 



B 



4.2 Double Born diagrams 

All 34 double Born graphs are presented in Figure [3j Graphs 1-7 are of topology 1, graphs 
8-16 are of topology 2, graphs 17-25 are of topology 3, and graphs 26-34 are of topology 4. A 
general diagram k with 1 < k < 34 equals 



M 



(2c) 



9 &1 £ 2 Xn,p 



d* 1± d* 2J .c fc r 1 fc l «/ 



(2c) 



J. 



(47) 



where C^ and rj-, 1 * 2 are identical to the single Born case, given in Eq. (45) and Eq. (46). But 
for double Born diagrams N min and N max are both larger by 1 in comparison to the single 



Born case. The maximum number of q— dependent denominators is still N n 



N min + 1 for 



topology 1 and A^ 



Amin + 2 for the remaining topologies. 



9 This similarity is very much expected due to the fact that Glauber exchanges do not change the large 
momentum fractions z 1 ,z 2 ;2 : 3 and, thus, the part of the amplitude that depends only on these fractions is 
identical to the vacuum case. Hence, the operators Ot k 
or 4 depending on the topology of the diagram. 



are equivalent to Eq. (18) and Eq. (19), and tk = 2,3 
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Figure 3: Double Born diagrams. See text for explanation of the different topologies 



The longitudinal integrals for double Born diagrams are defined similarly to Eq. (38) 
r (2c) _ f dgf dq 2 



2vr 2tt 



e i q ; 5zi+iq ; Sz2 Ag{Qu ~ i} ___ Ag (Q Nk: g Nh)j 



(4? 
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where qi is gi, q 2 or q 1 + q 2 dependent on the diagram. The results can be expressed through 



the same functions Ii, I 2 , h defined in Eq. (42)— Eq. (44) 



r(2c) 



(49) 



[Tiq^i Tlq 2 , ^g 12 J) 



h{tt 2 )/2, if rift = (1,0,1) 

I 2 (fi 2 ,ft 3 )/2, if n k = (1,0,2) 

J 3 (0 2 ,^ 3 ,^4)/2, if n fc = (l,0,3) 

12(^ + 02,^3), if n k = (1,1,1) 

/3(0 4 ,0 1 + 03,0 2 + 3 ), if rift = (2, 1,1) 
ia^fix + ^Ox + Qg), if rift = (1,2,1) 

13(^1 + ^2,^3,^4), if rife = (1,1, 2) 

with n gi being the number of q\ dependent denominators, etc. 



n gi + n q2 + n <?i2- F° r details on how to perform single and double Born longi- 



where n k = 

Obviously N t 

tudinal integrals see AppendixO The f2j in Eq. (49) are defined as: fix- • • f2 n<M for gi dependent 

propagators, £l nqi+ i. . . Q, nqi +Uq Jfor g 2 dependent propagators, and Vl 

for gi + q 2 dependent propagators. 

The color o pera tors e k for the double Born amplitudes in the basis of 24 basis elements 

J 2 ) 



n qi +n q2 +l- ■ ■ "n 9l +n g2 +n 9l2 



4.3) are provided in Appendix 



D 



^see section 

Even though we provided the rules necessary to evaluate any double Born graph in this 

B 



subsection, we also summarize all values for C k , U Pk , Pk , U Pfc , tPh , ± k 



I,„ in Appendix 



4.3 Squaring the matrix element 

In this subsection we combine vacuum, single and double Born amplitudes and derive a formula 
for the total squared matrix element averaged over the dense QCD matter. We start from 
general expressions for vacuum, single and double Born amplitudes 



M% = g 2 e\^ X , 



n,p?, 



M% = g 2 ^e 2 * X 



"■■PA 



A/f( 2c ) _ J ii is - 
Jvl n+2 ~ !J fc l fc 2 X.n,p 3 



J>f (a^ + c^E 3 ) J, 
3 

/"d^^e^^+^zS 3 ) J, 
J i 

J d*LL d$ 2± J2 e T ^ + 72? « 3 ) J- 



(50) 
(51) 
(52) 



The tensors structures ai j2 , /3i )2 , 7x j2 can be directly read off from the results in sections 2.2 , 4.1 



and 4.2 The basis of color operators for these three cases is: 



M 



,(i) 



(axa 2 )n, (a 2 ai) R , (53) 

(■■'■ = (aia 2 b) R , (ox6o 2 )fl, (6ox0 2 ) fl , (o 2 ax6)fl, (o 2 6ox)fl, (6o 2 Ox)a, (54) 

3 (2) = (aia 2 6i6 2 ) jR , (0x6x02^2)^, (Mio 2 & 2 )k, (a 2 aibib 2 )R, (a 2 &xax&2)-R> (&ia 2 aife 2 )ij, 

(aia 2 6 2 6i) R , (ai6i6 2 a 2 ) R , (Mx^g^k, (a 2 ai& 2 &i),R, (a 2 b 1 b 2 a 1 ) R , (6ia 2 6 2 ai) R , 
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{a 1 b 2 a 2 b 1 ) R , (a^b^R, {bib 2 aia 2 ) R , (a^ai&i)^ ( a 2&2Mi)fi, {b 1 b 2 a 2 ai) R , 
(b 2 aia 2 bi) R , (b 2 aibia 2 ) R , (b 2 bi<iia 2 ) R) {b 2 a 2 aib{) R , {b 2 a 2 biai) R , {b 2 bia 2 a{) R . 



(55) 



For brevity, we have omitted the overall medium color structure of (b)i for the single Born and 
{b\)i{b 2 )j for the double Born color basis elements. Combining, squaring, and averaging over 
the position of the medium scattering centers we get: 



spin, color 



,(0) 



f(l) 



,(2c) 



M™ 2 + M% 2 + M% 2 + ■ ■ 



(56) 



i± 



i 



g Tr -n-p 3 JJ 



p^ + w c TJ^^ Hq±)l2pl+Hq±W{ - q±)p ^ + ■■■]) 



( (o)t kA!) 



,( 2c )UA 2c ) 



Note that the term oc T^J^i n ^- 2 -Mn+2 vanishes and the term oc Tr_A/f„ +2 M-n+2 contributes 
to higher order in opacity [131 HI] . Here po is given by the vacuum splitting and has been 



calculated in section 2.3 The single and double Born terms p\ and p( 2c ) are: 

pi = E ( e i' } l e f) ( Re [Pis -fa + &VAu] J - Im lAVA* ~ Rffa] s3 ) > 

j'd 
P(2c) = 2 Y1 ( e f ] \ e ? } ) ( Re [ a i j' ' 7l J + a 2,j' -72 j] I + Im [a 2jjl -7i j - a*^ -72 J S 3 ) . (57) 



j >j 



In the equations above for the single Born expression we have used the fact that the Gram 
matrix of color basis vectors is symmetric. This is explicitly shown below in this subsection. 
The dot products between the tensor structures a,/3,7 indicate contractions, for example 
K-,-d\j = 77 ,■ ffi 18 !^?' The Gram matrices of the color vector basis necessary for the 
evaluation of the squared matrix element are straightforward to obtain 



^W 



T, 



R 



(efkf) 



T, 



R 



C\ c 2 c 3 c 2 c 3 c 4 

C 2 Ci C2 C3 C4 C3 

C3 C2 Ci C4 C3 C2 

C2 C 3 C 4 Ci C 2 C 3 

C3 c 4 C 3 C 2 Ci c 2 

C4 C3 C2 C3 C2 C\ 

C\ C 2 C3 C2 C3 C4 Ci Ci C2 C 2 C 2 C3 C2 Ci Ci C3 C2 C2 C3 C 2 C\ C4 C3 C2 
C 2 C3 C4 Ci C2 C3 C2 C2 C3 Ci C\ C 2 C3 C2 C2 C 2 C\ C\ C4 C3 C 2 C3 C2 Ci 



(58) 



(59) 

Because every element of these two matrices above is a number times a unit matrix in color 
space, the squared matrix element of both single and double Born amplitudes automatically 
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is a singlet in color space. The color factors as functions of the SU(3) Casimirs are 
a = C%, c 2 = C 2 F (C F - C A /2), c 3 = C F (C F - C A /2) 2 , 



C4 = Cf {C f - C A ){C F - C A /2) = 2c 3 - c 2 . 



(60) 



Unlike in vacuum, the squared two gluon amplitude in medium is in general not a singlet 
in Dirac space. This was also found for the single gluon probability kernel in Ref. [20]. In 
vacuum a\j and a^j are real and, hence, the E 3 piece cancels. This is not the casse in medium, 
because the longitudinal integrals have a non-zero complex phase. However, if the jet has been 
created by a pure QCD interaction, the trace Tr f | J J E 3 J = and the medium-induced two 
gluon emission factorizes from the production process, similarly to the single gluon emission 



4.4 Cascade approximation for the two gluon splitting function 

We define the 1 — > 2 and 1 — > 3 splittings in the presence of dense QCD matter 

2\ 2 9 2 



E 

spin, color 

E 

spin, color 



,(0) 



((1) 



(2c) 



AO + C + AC! + 



M 



(o) 

n+2 



,m 



m ' + Miy, + 



(2c) 



l n+2 



n+2 



q- 



q_ 



s 2 

s jk 



spin, color 



^(PwfaPk,*]) E \ M n\ 2 



>jkl 



spin, color 



The full splitting functions become a sum over the opacity series 

(p-^jk [pj ,Pk]) = (p!% k [pj ,Pk]) + (pQ jk [pj ,Pk\), 

(Pi^jki\Pj,Pk,Pi\) = (Pi%ki\PjiPk,Pi}) + ( p Qjki\PjiPk,Pi]}- 



(61) 

(62) 
(63) 

The first term corresponds to the vacuum splitting functional and the second term corre- 
sponds to the first order in opacity term, including both single and double Born graphs. Note 
that Pi,Pj,Pk,Pi are momenta of external partons and independent of the medium averaging, 
which is not shown in the above equations but is present in the second terms. The medium- 
induced cascade formula, similarly to the vacuum case, is based on the approximation that the 
probability to emit two gluons can be approximated by a product of single gluon emissions. 
This approximation is valid up to certain interference terms. The splitting function of the 
medium-induced "cascade" is 



/ pease 



\PUP2,P3}} [ ' = 



Sl23 i 



[Pi% q \p 2 ,Pl +P3])(Pi% q \pi,P3}) + (Pi% q [P2,Pl +Ps]){Pi% q [pl,P3}) 
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10 It is identical to it due to our normalization 
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S23 
(P^ OT [pi + P2 , P 3 ] ) (P<S<7<7 [Pi , P2] ) + (Pq\q [Pi + P2 , P 3 ] ) (-PJ2 M [Pi , P2] ) , , , 

+ j . (64) 



This equation is derived analogously to the one in vacuum, Eq. (32), and takes into account 
that the interaction with the medium can happen either in the first or the second splitting;] 
The medium-modified 1 — > 2 splitting functions are related to the medium-induced splitting 
kernels x dN/dx dkj_ calculated in Ref . [201 EI] arid are reviewed in Appendix [A] 

5 Angular distributions of splitting functions 

In this section we study the angular distributions of the collinear vacuum and medium-induced 
splittings. We start with an overview of coherent branching and angular ordering, following 
closely Ref. [15]. Consider an arbitrary hard process with a total of n incoming and outgoing 
quarks and/or gluons and an exclusive differential cross section a n . In addition, we define 
cr n+ i as the lowest order differential cross section to emit an ultrasoft (eikonal) gluon with 
momentum scaling (A 2 , A 2 , A 2 ) from either of the external legs. Using the well known eikonal 
approximation of QCD we find 

da n+ \ = da n > C;~ Wa, (65) 

+ u 2vr 2vr ^ 3 v v ; 

where u is the energy of the emitted gluon, Cy is a color factor and 

w " 2 p*-Pi = 1 - cos ^ . (66) 

1 Pi-<lPj-<l (1 - COS iq ) (1 - COS Ojq) 



Each term of the sum in Eq. (65) corresponds to a different interference term, where the 
ultrasoft gluon is attached to the leg i in the matrix element and leg j in the complex conjugate 
of the matrix element. Thus, the ultrasoft branching depends on the global structure of the 
event. We have assumed that all external legs are massless, pi is the momentum of leg i, q is 
the momentum of the emitted gluon and the angles between legs i and j and leg i and the soft 
gluon are defined as 0^ and 8 iq , respectively. The function Wij has the well known property 
of angular ordering. Namely, if one rewrites 



where 



W V = W$ + W®, (67) 



W $ = \ i W *i + i 1 a ~ \ ^—r] ' ( K 

J 2 V 1 — cos fag 1 — cos "jq J 



11 There can be medium interactions in both splittings at higher orders in opacity. 
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and a similar definition for Wf- with i •H- j, then W/ has the property 



i.i 



d(f>i q 
2tt " ij ~ 1- cos ft- Wl "" " ;,/ 



WS- , .__„ Q(%~^)- (69) 



The integration over the azimuthal part of dQ in Eq. (65) is performed while fixing the z 
axis along the direction of parton i. Thus, the angle 8 iq is kept fixed while 9j q varies as 
cos 9 j g = cos 6 ij cos 6 iq + sin 8ij sin 9 iq cos <fii q . Eq. (69) means that after azimuthal averaging 



the interference term W\- emits radiation only inside the cone R\\ with opening angle 6{ 



centered around parton i. W\f obeys an equation analogous to Eq. (69) and only radiates soft 



gluons inside the cone R^ . In the remainder of this paper we will refer to the cones R^ and 



>[./] 



Rl as angular ordered cones. It follows from these properties and equation Eq. (65) that in 
the eikonal approximation the radiation obeys angular ordering in the sense that the emitted 
gluons are emitted only inside the angular ordered cones for all i and j. 

Having reviewed the known properties of coherent branching, we turn to the question of 
how the collinear branchings behave in terms of angular distributions. In the high energy 
factorization picture of hard scattering processes there are three widely separated distance 
scales: the scale of the hard process ~ l/\/s, the scale of collinear splittings or showering ~ 
1/£>_L> where p± is the scale of the transverse size of the jet, and the scale of soft recombination 
processes ~ 1/Aqcd- The coherent branchings discussed in the last paragraph which lead 
to angular ordering correspond to the soft scale. Angular ordering is widely used in the 
literature on parton showers and is implemented in some of them, for example HERWIG 
and PYTHIA. To our knowledge, the study of angular distributions of collinear splittings, 
which are characterized by the intermediate scale of the parton shower, does not exist in the 
literature. In the following we perform such a study in vacuum and medium using full 1 — > 3 
collinear splittings which include all interferences. For the vacuum we use results from Ref. 
[8] , where all 1 — > 3 parton collinear splittings have been calculated. In section [2] we derived 
one of these splittings, q — > ggq, in SCET and confirmed the result in Ref. [8]. In medium we 
use our new SCETq result for the q — > ggq splitting presented in section [4] 

For later use, we define 

Xij = ~\ (Wij - l — r - 1 -—) , (70) 

2 V 1 — COS "iq 1 — COS Vjq J 

which has the property of anti- angular ordering 



d(pi Q y _ 1 

2tt lJ 1-cosft wv "'" "" 



X U = ~, 7T ©(^ - %)■ ( 71 ) 



'iq 



5.1 The vacuum case 

We consider the five 1 — > 3 splitting functions, calculated in Ref. [S], and study their an- 
gular distributions. Our goal is to clarify if collinear splittings exhibit a feature like angular 
ordering of soft coherent branching discussed in the previous subsection. It is clear that the 
notion of angular ordering is only applicable in deterministic parton showers with sequential 
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branching. The obvious example is precisely the q —¥ ggq splitting, where the two gluons are 
indistinguishable. In order to define a notion of first and second splitting, we choose the limit 
when one of the three partons in the final state is much softer than the two others, z\ <C z 2 ,z 3 . 
One has to be careful in taking this limit, since the collinear power counting breaks down if 
parton "1" becomes too soft. So we have to ensure that the energy carried by parton "1" is 
much smaller than the one carried by partons "2" and "3", but still much larger than Aqcd- 
Or, in other words, we take the limit in the collinear branching 1 — > 3 such that the second 
branching is at larger distance than the first one, but still in the collinear region, not in the 
soft recombination regime. In contrast to the ultrasoft branching in the previous subsection, 
this limit is process independent because the collinear splitting functions we started with are 
process independent. 

Taking this limit for the five splittings: q — > q'q'q, q —¥ qqq, q —¥ ggq, g — V gqq, g — > ggg 
yields 



V <?<)->9i9 2 93/ 
\-^9o-»9l9293/ 



C F T R (1 - c 23 ) (1 - z 2 ) (2(1 - z 2 ) + zl 



CfQ. - c 23 ) 

Zl 



Z-2 



w w + x 23) . 



(72) 



T.( (1 -^ (1 + (1 -^ 2) (^ + X 2 3 



V 90 -*-91 92 93/ 
V. 90->-9l9293/ 
'-^90-^919293/ 



4CW1 - c 23 ) 



zl 




2T R (1 - 


C23) 


zl 




4Ci(l - 


-c 23 ) 



z 2 (l- z 2 



+ "-f^ (Vg + X 23 )) + 2(C F - C A /2) X, 
I-Z2 + z\j1 



Z-2 



C F (W^+X 23 )+C A (W, 



23 1 



22(1 - z 2 ){zi + (1 - z 2 ) 2 ) C F W% + W£ + C A (X- 



M 



r [3] 



^23, 



^2(1 - z 2 \ 



z 2 



1-Z2 



- — - + z 2 (\ - z 2 )\{wtl + Wf 3 ] - X-r, 



We see from the presence of both terms W\ l j and X^ that the splittings are neither angular 
ordered nor anti-angular ordered, however, some individual pieces are. For example, the non- 
abelian part of the third splitting and the abelian part of the fourth splitting are angular 
ordered. The identical particle piece of the second splitting proportional to C F {C F — Ca/2) 
and the non-abelian part of the fourth splitting are anti-ordered. All other pieces are neither 
ordered nor anti-ordered. Note that, as expected, the last three splittings are proportional to 
the reduced vacuum 1 — > 2 splitting of the initial parton into partons " 2" and " 3" . Conversely, 
the first two splittings do not exhibit any similar relation. In addition, the soft behavior of 
the first two splittings compared to the last three differ. The first two splittings, where the 
anti-quark is taken to be the softer parton, are proportional to 1/u ~ l/^i and thus free of 
soft singularities after including the phase space factor dz\Zi. The last three splittings are 
proportional to 1/zf which leads to a soft singularity. 
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Figure 4: Angular distributions of 1 — > 3 splittings in the vacuum. 
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Next we compare the small z\ limit of the full q —¥ ggq splitting, given in Eq. (72), to the 



small Z\ limit of the cascade Eq. (32) 



V*VnnW - U - C23J— ^^2(1 ~ Z 2 ) 

Z\ z 2 

x (C F (W| + X 23 ) + C A (V| ] + X 23 ) ) . (73) 

The cascade reproduces the abelian part exactly, while only the singular behavior is reproduced 
for the non-abelian part. The numerical comparison of the full splitting to the cascade is shown 
in the middle right panel of Figure |4| As expected from the analytical formulas, the cascade 
reproduces both collinear singularities. In the tail of the distribution the cascade is larger by 
a factor (Cp + Ca)/Cf compared to the equivalent piece (X23) in the full splitting. 

To visualize the angular distribution we plot the splitting function averaged over the az- 
imuthal angle of the softer parton (z\) with respect to the decaying parton as a function of the 
angle between the decaying parton and the softer parton (z\). No other phase space factors 
are included. We fix the angle between the second (third) and the decaying parton to be 10 
degrees (20 degrees). We set the energy of the initial quark E = 100 GeV, z\ = 0.03 and 
z 2 = 2/3, Z3 = 1/3. Note that z\ is much smaller than z 2 and Z3 but the enrgy of parton "1" is 
still much larger than Aqcd- Thus, this choice of parameters obeys the desired limit. For the 
numerics we exploit that partons "2" and "3" have to be approximately back-to-back in the 
transverse plane in order to balance the total transvers momentum. In Figure [4] we present 
these plot for all five splittings, which we discussed in the following: 

• Each splitting has a collinear singularity at angles of 10 and 20 degrees as expected. 

• The steepness of the angular distributions outside the angular ordered cones (at 40 and 
50 degrees with respect to the initial decaying parton) is ~ 1/^og i n the cases when there 
is angular ordering and ~ 1/^oq when there is no ordering, consistent with equation 



Eq. (72). Note that in the situations when there is ordering, radiation is still present 
outside the angular ordered cones, though it is power suppressed compared to the case 
of no ordering. The reason for this is that we are averaging over the azimuthal angle 
with respect to the initial (decaying) parton, not with respect to one of the partons "2" 
or "3" like in Eq. (69). Since W}j and W\y become non-positive definite outside the 



angular ordered cones they cannot be interpreted as probabilities anymore. Since the 
solid angle of the emitted gluon is ~ d^, the amount of radiation outside of the cone 
depends on the cone size logarithmically for the power law ~ l/^og> but by an inverse 
power for the power law ~ 1/6, 



o q - 



In the middle left panel of Figure [4], for the go - > 9i92Q3 splitting, the abilian contribu- 
tion dominates in the tail, while the non-abilian contribution to the full result is only 
marginal. This can be understood qualitatively. When a gluon is emitted at large angle, 
it cannot resolve the small angle between g 2 and q^. Thus, the gluon is effectively emit- 
ted from an on-shell quark (q ). A similar argument holds for the splitting g — > g\q 2 qz 
and leads to the non-abilian radiation dominating over the abilian radiation in the tail, 
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consistent with the bottom left panel of Figure |4j The same qualitative analysis yields 
that the tail of the angular distribution for the splitting g — y ggg is not ordered. 

• Note that we plot the absolute value of the angular distribution for the identical-particles 
term of the splitting q — y qqq and the abelian term of the splitting g — y gqq. In both 
cases, the distribution is not always positive and could not be shown in a logarithmic plot. 
In the first case the true contribution is positive for #01 < 10°, zero for 10° < #01 < 20° 
and negative for #01 > 20°. In the second case, the true contribution is negative for 
#01 < 10°, zero for 10° < 9 01 < 20° and positive for 9 01 > 20°. The difference in sign 
between first and second case is due to Cp — Ca/2 < 0. 

In conclusion, we find that inside the collinear parton shower there is no angular ordering 
in contrast to the ultrasoft coherent branching. The implications of this result for parton 
showers remain to be studied phenomenologically by concentrating on observables related to 
collimated, isolated jets. This will be done elsewhere. 

5.2 The dense QCD matter case 

In section [4] we calculated the medium-induced splitting q — y g±g2q3 using SCET G , which 
was the most technically demanding part of this paper. See Appendix [X] for details on the full 
calculation and an approximate reduced formula valid in the small Z\ limit. In this section, 
we perform an analysis of the angular distributions in the medium-induced splitting similarly 
to the one in vacuum. In vacuum, we considered five splittings and were able to analize 
the angular distributions in the small Z\ limit analytically. In medium, we just perform a 
numerical analysis of the splitting q — y ggq. 

To model the QCD medium, we use the following input parameters: the Debye screening 
scale in the medium is fi = 0.75 GeV, the size of the medium L = 5fm, and the elastic 
scattering length of gluons in medium X g = 1 fm. These values have been used in [21] and are 
characteristic of the quark-gluon plasmas created at RHIC and LHC. The numerical results are 
shown in Figures [5j [6j [7j These plots are the medium equivalent of Figure |4| but this time we 
consider two different sets of parameters. In scenario 1 (Figure [7] and top plots in Figures [5j |6| 
we use the similar values like in section |5.1| for vacuum: E = 100 GeV, z\ = 0.03, z 2 = 



0.643, #20 = 10°, #30 = 20°. In scenario 2 (bottom plot in Figures § [§ we use: E = 100 GeV, 
z 1 = 0.03, z 2 = 0.282, # 20 = 25°,# 30 = 10°. We present the total medium splitting (solid 
black curve), medium cascade (dot-dashed green curve) and the vacuum splitting (dashed red 
curve). In Figure [5] we compare medium to vacuum splitting, in Figure |6J medium splitting 
to the medium cascade and in Figure [7] medium to vacuum splitting but in three dimensions. 
We make the following observations: 

• The collinear singularities are present in both single and double Born graphs. When 
combining single and double Born graphs we find large (90% to 99%) cancellation for 
both scenarios. The collinear behavior, corresponding to the gluon z\ being parallel to 
the quark z 3 , is in both scenarios significantly reduced. 

• For both scenarios the tail of the angular distribution is larger for the medium-induced 
splitting than for the vacuum splitting. Moreover, in the direction of the parent parton 
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(small 6*01) there is significant cancellation in the splitting probability. These features are 
in agreement with the previously noted features of medium-induced radiation. Namely, 
the 0(a s ) q — >■ qg splitting in dense QCD matter is wider than in the vacuum p2 



As one can see in Figure |6j the cascade formula for the medium-induced splitting de- 
scribes the qualitative features of the full splittings for angles between the peaks and in 
the tail reasonably well. The cancellation of the splitting probability along the direction 
of the parent quark is not reproduced by the cascade. 

We compared our full medium splitting formula to the approximate formula presented 



in Appendix A. 2 which includes only topologies two and four and is valid for small Z\. 
For both scenarios the difference between the exact and approximate formula is smaller 
than the visible thickness of the lines in our logarithmic plot, Figure [5j This is a nice 
cross check on our numerics. 

• In Figure [7] the spacial distribution of the medium- induced splitting is compared to that 
of vacuum. The collinear radiation is significantly reduced in the medium, while the 
medium-induced radiation is larger in the far tail. 

In conclusion, the medium-induced splitting exhibits no angular ordering or angular anti- 
ordering, similarly to the vacuum splitting. The splitting probability distribution is larger in 
the tail in comparison to the vacuum splitting. Moreover, there is a cancellation of this prob- 
ability in the direction of the original parton. These features have been described previously 
for the lowest order medium-induced parton branchings. From a practical point of view, our 
results imply that in constructing Monte Carlo generators to describe jet physics in heavy ion 
collisions, an approach without angular ordering would be preferred. 
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Figure 5: Angular distributions of 1 — > 3 splittings in vacuum and the medium. Color coding 
is as follows: red— vacuum splitting, black— medium splitting. The first plot corresponds to 
scenario 1, the second plot to scenario 2. Further details are given in the text. 
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Figure 6: Angular distributions of 1 — > 3 splittings in the medium. Color coding is as follows: 
black— medium splitting, green— medium cascade. The first plot corresponds to scenario 1, 
the second plot to scenario 2. Further details are given in the text. 
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Figure 7: Angular distributions of 1 — Y 3 splittings in vacuum and medium. Further details 
are given in the text. 

6 Conclusions 

In this paper we studied the the final-state angular distributions of 1 — > 3 collinear splitting 
functions at order 0(a 2 s ) in vacuum and in dense QCD matter. We concentrated on the split- 
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ting q — > ggq. By comparing the angular distribution of this splitting function to the ones of 
the other splitting functions, for instance g — > ggg, we showed that q — > ggq is a representa- 
tive example. Using SCET, we calculated the vacuum splitting function to demonstrate that 
collinear modes are sufficient to reproduce the result obtained by Catani and Grazzini [8]. In 
dense QCD matter we used SCETg to derived the medium-induced q — ¥ ggq splitting to first 
order in opacity, keeping the full Z\, z 2 dependence. 

In vacuum we studied all five 0(a 2 s ) splitting functions, originally derived in Ref. [8]. In all 
cases we find no angular ordering. This result is also supported by the fact that the angular 
distribution of the q — > ggq splitting can be reproduced well by a parton cascade based on 
binary branchings. Note that our approach differs from the coherent branching ansatz which 
yields angular ordering and is widely used in parton shower generators like HERWIG. We 
study angular distributions of collinear splittings instead of ultrasoft gluon emissions. Thus, 
the momentum scaling in our approach differs from the one in coherent branching and hence, 
the different result is not unexpected. Nevertheless, the qualitative argument of color screening 
used to explain Sudakov effect in QED still works in our case. In large angle gluon emission 
from a qg antenna, the emitted gluon is only sensitive to the color charge of the initial quark. 

Our results indicate that the proper angular distributions inside the collinear parton shower 
are different from the traditional coherent branching ansatz. Nevertheless, the traditional ap- 
proach to parton showers, which applies angular ordering to the collinear splitting functions, 
is claimed to resum large infrared Sudakov logarithms and is phenomenologically successful. 
One thing which would be interesting to verify is whether the amount of collinear radiation 
leaking outside of the angular ordered cones leads to a significant correction to parton shower 
phenomenology. This would have to be checked for intra-jet observables, such as jet shapes of 
well-isolated jets. We leave this for future work. The same conclusion holds for the medium 
induced parton shower. Our detailed analysis found no evidence of angular ordering or an- 
gular anti-ordering. An important feature is that the noticeably broader angular distribution 
compared to vacuum, found in 0(a s ) 1 — > 2 branchings, persists to higher order. 

A Medium- induced splitting functions 

In this appendix we review the basic formulas for medium-induced splitting functions and how 
they are related to the splitting kernels. 

A.l Leading order 0(a s ) splittings 

The medium-induced splitting kernels have been calculated in Ref. J20J EI] retaining the full 
x dependence (beyond the soft gluon approximation) using SCET G . We rewrite them in a 
slightly more compact form: 

diV"W a s „/,•-,, . f dAz l9 1 dcr c i v^ ((\,« ^ \ /„.% 

-r-^rr = 77^P^c(x) / YTT d C ^ — ^E a ° l-cos$ fc . 74 

dxd 2 k_L 2tH J Xi(z) creid 2 q_L^ 

We define the transverse vectors 

A_L=k_L, B ± = k ± +a;q ± , C ± = k_L - (1 - x)q ± , D ± = k ± - q ± . (75) 
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In terms of these vectors, the five phases $^ are equal to 

$1 = S>B 2 ± , $ 2 = VC 2 ± , $ 3 = *(C^ - B 2 ± ), $ 4 = VA 2 ± , $ £ 

A* 
where \l/ 



^(A 2 



Di), 



x(l — x) n-p 



(76) 



We recall that Az is the spacial separation between the hard scattering, producing the collinear 
parent parton, and one of the subsequent medium interactions via a Glauber gluon exchange. 
The coefficients at are summarized in the following table: 



k 


a ( r q9) 




a [ r m 


1 


b-(b-c+^) 


2b(b-a- s=5) 


2b.(b-a+^) 


2 


c(2c-a-b) 


2c-(c-a-^) 


2c- (c a+^) 


3 


be 


be 


r> b-C 


4 


a(d — a) 


a(d — a) 


a(a— d) 


5 


—ad 


—ad 


r> ad 

Z 7V r 2 -l 



where a = A^/A^ and b, c, d are defined similarly. To relate the medium-induced splitting 
function to the splitting kernel we have to remove the phase space contributions; recall that 



dJVW 
dcr n+ i = d<7„ - — — — dx d 2 k i 



d$ 



n+l 



d$ r 



dx d 2 kj 
1 



x(l — x) IGtt 3 



dx d k i . 



(77) 
(78) 



Using that the splitting function is defined as ratio between the squared matrix element after 
emission and before emission, gives 



(P {1) )il 



)W, 



dAz 2 1 dad v^ u). . 



PW T b 

Ai(z) a e id 2 q 



(79) 



fe=i 



Note that the q — > gq kernel is obtained from the q -^ qg kernel via the substitution x — > 1 — x. 
For the normalized elastic scattering cross-section we get 



1 da, 



cl 



l>- 



Cei d 2 q± TV (q 2 ± + /i 2 



,2' 



(80) 



consistent with section HI 



A. 2 Next-to-leading order 0{a 2 s ) splittings 

For the q — > ggq splitting function to first order in opacity we find 



(P W ) 



q^-ggq 



1 s 2 

1 A 12 3 

N c 4 



-z 3 



dAz 
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d 2 qj . 



1 da, 



cl 



a ei d 2 qj 



(Pi + P(2c)) 



where p\ and p( 2c \ are provided in general form in Eq. (57). For the special case when parton 



one has softer momentum than the other two partons, z\ <C z 2 , z 3 , these formula reduces to 

* « 4 (i -* 2 + 1) £ <4I , '|r«"l4 1, ')^;>c« (uf > -uf >) (uf > .uf > 



2 , 

k'k=l 



xRe /<;'•/«, 
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(2 \ lu 

i - « + f ) E <4^ir« 2 'l4 2) ') cf cf (uf > -uf ») (uf ' -uf » 
' fc'=l,2;fc=l,18 

x2Re4 2c) , (82) 

where for pi the sum over k runs over the 10 single Born graphs of topologies 2 and 4; for 
P(2 C ) the sum over /c runs over the 18 graphs of topologies 2 and 4; as well in vacuum only 
topologies 2 and 4 matter in this limit. Our notation for U*- 11 - 1 , U*- 12 - 1 , U*- 21 **, U < - 22 ' ) is as follows. 
For a given single Born graph, \j( u \\j( 12 ) are the first and second transverse vectors in the 
corresponding entry in the third raw in the Table |l| Similarly, \j( 21 \\j( 22 ') are the first and 
second transverse vectors in the corresponding entry in the third raw in the Table |2J T^ 
and r^ 2 ) are defined in Eq. (58) and Eq. (59). All other ingredients follow from the rules in 



section [4] for longitudinal integrals, coefficients C^ and color operators e^. 

B Feynman Graphs 

In section [4] we presented rules to extract an analytic expression for any single and double Born 



Feynman diagram. Eq. (45) and Eq. (46) are valid for any Feynman graph. In this appendix 



we give explicit expressions for all coefficients C&, longitudinal integrals Ik and transverse 
vectors U Pfci , p ,U P Pfc4 . 

In Figure [2] and [3] all 19 single Born and 34 double born Diagrams are shown, respectively. 
The corresponding analytic components can be read off Table [T] and [2] The frequencies are 
defined as 

Q = Q(p! + p 2 +]?3,q±), 

fii = £l(pi, qjj, Q 2 = Q(j>2, q±) ^3 = 0(p 3 , qj_), 

tt 4 = tt(p 2 +p 3 ,q ± ), n 5 = n(p 1 +p 3 ,q ± ) tt 6 = tt(p 1 +p 2 ,q ± ), (83) 

and 

ft, = fi,(<u -)• 0), Qi = fi,(<u -► -qj, where Z = 0, 1, 2, 3, 4, 5, 6. (84) 
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Table 1: Entries for single Born graphs. 



29 



k 


c k 


XJii \jn 
Pk 1 ,Pk 2 Pk S ,Pk 4 


/r j /(-^) 


1 


1/«123 


— 


/ 1 (n )/2 


2 


l/n-p 


— 


/ 2 (n ,n 2 )/2 


3 


1/n-po 


— 


7 2 (fi ,^i)/2 


4 


1/n-po 


— 


7 2 (O ,O 3 )/2 


5 


1/n-po 


— 


/ 2 (Oo,0 2 + 3 ) 


6 


1/n-po 


— 


/ 2 (Oo,Oi + 3 ) 


7 


1/n-po 


— 


/ 2 (Oo,0 2 + Oi) 


8 


1/Sl3«123 


T P 1 T P 2 

Pl,P3 KJ P2,P3+Pl 


/i(n )/2 


9 


l/n-p si3 


T P 1 T P 2 

Pl,P3 P2,P3+Pl 


I 2 (O ,O 2 )/2 


10 


±/(n-p ) 2 ( Zl + z 3 ) 


T P 1 T P 2 

P1:P3 P2,P3+Pl 


/ 3 (Qo,^i,^ 5 )/2 


11 


±/(n-p ) 2 ( Zl + z 3 ) 


T P 1 T P 2 

P1,P3 P2,P3+Pl 


I 3 (O ,O 3 , 5 )/2 


12 


l/n-p s 13 


T P 1 T P 2 

^Pl,P3 ^P2,P3+Pl 


/ 2 (O ,n 5 )/2 


13 


^-/(n-p ) 2 (z 1 + z 3 ) 


Tpl Tp2 

pi,P3-q P2+q,P3-q+pi 


73(00,^2 + ^3,^2 + ^5) 


14 


l/(n-Po) 2 (zi + z 3 ) 


TjJl TJJ2 

Pl-<?>P3+<? P2,P3+Pl 


/ 3 (o , o 5 , o : + o 3 ) 


15 


l/(n-p ) 2 (^i + z 3 ) 


TjJl TJJ2 

^Pi-g.pa v -'p2+<7,P3+Pi-<? 


7 3 (O , 2 + Oi,0 2 + 5 ) 


16 


l/n-p s 13 


T P 1 T I- 72 

Pl,P3 ^P2+(?,P3+Pl-l3 


/ 2 (O ,O 5 + O 2 ) 


17 


l/s 23 s 123 


T P 1 T I- 72 

P2,P3 ^Pl,P3+P2 


/i(O )/2 


18 


l/n-p S23 


T P 1 T P 2 

P2,P3 Pl:P3+P2 


/ 2 (n ,no/2 


19 


^/(n-p ) 2 (z 2 + z 3 ) 


T P 1 T P 2 

P2,P3 Pl,P3+P2 


/3(Oo,o 2 ,o 4 )/2 


20 


l/(n-p ) 2 (z 2 + z 3 ) 


T P 1 T P 2 

P2,P3 Pl,P3+P2 


/3(no,n3,o 4 )/2 


21 


±/n-p s 23 


T P 1 T P 2 

P2,P3 Pl,P3+P2 


J 2 (O ,O 4 )/2 


22 


^/{n-p ) 2 (z 2 + z 3 ) 


Tpl Tp2 

V ^P2,P3-9 V ^Pl+g,P3-g+P2 


7 3 (O , Oi + 4 , Oi + 3 ) 


23 


^/{n-p ) 2 (z 2 + z 3 ) 


IP 1 IP 2 

P2-<?,P3+g^Pl,P3+P2 


J 3 (Oo,04,02 + 03) 


24 


l/(n-p ) 2 (z 2 + z 3 ) 


IP 1 IP 2 

^P2-<?,P3 V ^Pl+9,P3+P2-(? 


/3(Oo,Oi + 02,Oi + 4 ) 


25 


l/n-p S23 


Tpl TP2 

P2,P3 V ^Pl+9,P3+P2-'/ 


/ 2 (Oo,0 4 + Oi) 


26 


I/S12S123 


T P 1 T P 2 

Pl,P2 ^Pl+P2,P3 


/i(O )/2 


27 


l/(n-p Q ) 2 (zi + z 2 ) 


T P 1 T P 2 

Pl,P2^Pl+P2,P3 


/ 3 (Oo,0 2 ,0 6 )/2 


28 


V(^-Po) 2 (^i + z 2 ) 


T P 1 T P 2 

P1,P2 P1+P2.P3 


7,(00,0!, 6 )/2 


29 


l/n-p s 12 


T P 1 T P 2 

P1,P2 P1+P2.P3 


/ 2 (Oo,0 3 )/2 


30 


l/n-p s 12 


T P 1 T P 2 

Pl,P2 P1+P2.P3 


J 2 (0 ,0 6 )/2 


31 


l/{n-p ) 2 (z 1 + z 2 ) 


Tpl Tp2 

Pi,P2-q ^ pi+P2~q,P3+q 


/ 3 (O ,O 3 + O 2 ,O 3 + O 6 ) 


32 


l/(n-p ) 2 (^i + z 2 ) 


IP 1 IP 2 
pi-q,P2 ^pi-q+P2,P3+q 


/ 3 (Oo,o 3 + n 1 ,n 3 + n 6 ) 


33 


l/(n-p ) 2 (zi + z 2 ) 


IP 1 IP 2 

Ky pi~q,P2+q Ky Pl +P2 ,P3 


7 3 (O , 6 , Oi + 2 ) 


34 


l/n-p Q si 2 


T P 1 T P 2 

Pl,P2 Pl+P2-g,P3+9 


/ 2 (Oo,0 3 + 6 ) 



Table 2: Entries for double Born graphs. 
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C Longitudinal integrals 

In this appendix we derive all necessary formulas to calculate the single and double Born 
longitudinal integrals that appear in our paper. All single Born diagrams are of the form 



dq 



Hi, 



/,?>,, fcH/^-n^, (85) 

where n q is the total number of g-dependent propagators in the graph and «» = Q(Qi, q_i_) + 
ie/n-Qi. Integration is straightforward using Cauchy's theorem 

n q n q 

lg q \a i ,6z>0) = (-i)Y i e(1m ai )e ia * s ' J] — — , (86) 

i=l l=l,tyi l % 

/«(«,, Sz = 0) = 5 lim o lW( ai , Sz), if n q > 1, (87) 

j£>(«i, 8z = 0) = ± 5 ]im o /W(«„ Sz), if n q = 1. (88) 

The first two equations simply follow from Cauchy's theorem when closing the contour above. 
The third case is more subtle since the boundary term at infinity cannot be neglected which 
causes the factor 1/2. We explain this in more detail at the end of this section. The single 
Born diagrams in this paper take values of n q — 1, 2, 3 only. The three corresponding master 
formulas for these longitudinal integrals are 

I^ica, Sz) = -ie taiSz , (89) 

„ia 2 5z _ iaxSz 

l£\a u a 2 , Sz) = i , (90) 

II \a 1 ,a 2 ,a 3 ,Sz) = i I . (91) 

\ a 2 — ai a 3 — «i / «3 — «2 

Any double Born integral can be written as 

(n qi \ / n ?2 

/ n 912 , \ 

x II ^ » C 92 ) 

where a^fij,^ are frequencies that appear in the poles of A g (Qi,q) with Glauber gluon 
momenta q±,q 2 , q\ + q 2 , respectively. Performing the q^ integration using Cauchy's theorem, 
the remaining q[ integration can be expressed by a single Born integral 

l ni ,„„,,„„ («<; &5 7fc, <&*i, ^2) = (93) 



'91 '"92' 93 
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fe=i 
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lm ~lk 



J n„ +n q , (Oi! 7fc - ft, fol - fofe). 



This equation solves any longitudinal integral of the double Born graphs in this paper. 
For the first order in opacity calculation we need the contact limit of this integral, i.e. 



r(2) 



Eq. (93). 



(a*; /3j] 7fc, 8z, 5z). The result can be obtained by applying Eq. (86)-Eq. (88) to 



The longitudinal integral of any double Born diagram in this paper can be relateed to one 



of seven master integrals by n k 



\Tlq 11 Tlq 2 , Tlq 12 



I. This master integrals can be calculated to 



r(2c) 



-i)-< 



Ji( 7 i)/2, if n k = (1,0,1) 

/ 2 ( 7l , 72 )/2, if n k = (1,0,2) 

4(7i,72, 73) A if n fc = (l,0,3) 

I 2 {ocx + /3i,7i), if rife = (1, 1, 1) 

^(ai + /3i,a 2 + /3i,7i), if n fe = (2, 1, 1) 

h{ai + Pi,ai + /3 2 ,7i), if ™fc = (1,2,1) 

J 3 (ati + /3i, 71,72), if n k = (1,1,2) 



> ■ 



(94) 



where the functions ii, I2, I3 are the same as in Eq. (89)-Eq. (91 ). There is one technical issue 



when deriving Eq. (94). By definition at,/3j,7fc all have positive imaginary parts, due to the 
ie-prescription. However, one has to be careful when considering the complex number 7^ — (3j 



in Eq. (93). For all calculation we are concerned with this number has a negative imaginary 



part. Indeed, we have checked that the results of the double Born integrals do not depend on 
the sign of the imaginary part of 7^ — 0j. The crucial point is that none of these numbers lie 
on the real axis. 

Now let us return to the subtlety in I{ (an., Sz — 0). For n q = 1 and 8z = the boundary 
term at infinity does not vanish and the integral /} («i, 5z = 0) becomes ill-defined. To solve 
this issue one has to recall that the factor v(q) also depends on q~, v(q~,q±) ~ 2 2 l , - )2 , 

but we dropped this dependence since it is power-suppressed. Re-introducing this extra q~ 
dependence cures the boundary term probem and yields l4, {on, 5z = 0) = —i/2 and the overall 
factor v(0, q_i_). One can convince oneself that in the other integrals, like /„ >i(<2j, 8z = 0) or 

In q (on, 8z > 0), re-introducing the factor v(q~, q_i_) does not change anything to leading power 
in the EFT. 

D Color operators 

In this appendix we present the color operators of all single and double Born graphs in the 
paper. For the single Born graphs the color operators in terms of the basis vectors given in 
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Eq. (54) are 
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For the double Born graphs the color operators in terms of the basis given in Eq. (55) are 
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